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The aim of this paper is to derive the averaged governing equations for non-degenerated 
oscillatory flows, in which the magnitudes of mean velocity and oscillating velocity are 
similar. We derive the averaged equations for a scalar passive admixture, for a vectorial 
passive admixture (magnetic field in kinematic MHD), and for vortex dynamics. The 
small parameter of our asymptotic theory is the inverse dimensionless frequency l/cr. Our 
mathematical approach combines the two-timing method, distinguished limits, and the 
use of commutators to simplify calculations. This approach produces recurrent equations 
for both the averaged and oscillating parts of unknown fields. We do not use any physical 
or mathematical assumptions (except the most common ones) and present calculations for 
the first three (zeroth, first, and second) successive approximations. In all our examples 
the averaged equations exhibit the universal structure: the Reynolds-stress-type terms 
(or the cross-correlations) are transformed into drift velocities, pseudo-diffusion, and two 
other terms reminiscent of Moffatt's mean-fields in turbulence. In particular, the averaged 
motion of a passive scalar admixture is described only by a drift and pseudo-diffusion. 
The averaged equations for a passive vectorial admixture and for vortex dynamics possess 
two mean-field terms, additional to pseudo-diffusion. It is remarkable that all mean-field 
terms (including pseudo-diffusion) are expressed by invariant operators (Lie-derivatives) 
which measure the deviation of some tensors from their 'frozcn-in' values. Some physical 
assumptions and the results can be build upon obtained averaged equations. Our physical 
interpretation suggests purely kinematic nature of pseudo-diffusion. 



1. Introduction 

Oscillatory flows play important roles in all branches of fluid dynamics, especially in 
wave theory, in biological and medical fluid mechanics, in geophysical and astrophysi- 
cal fluid dynamics. If a chosen fluid flow is not oscillating, it can be easily made so by 
adding some oscillations to boundary conditions, etc. There is no unique definition of an 
oscillatory flow. In this paper we suggest that such a flow contains velocity oscillations 
with frequency higher than inverse characteristic times of all other co-existing motions. 
It can be expressed as a ^ 1, where a is dimensionless frequency of oscillations. Our 
term a non- degenerated oscillatory flow means that the mean and oscillatory velocity 
are similar in magnitudes. This class of flows is rather restrictive: in many applica- 



tions the oscillatim 


I part of velocity is dominating (hence, the flows are degenerated), 
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Stokes (1847)[ 


Craik (1985) Magar & Pedley (2005)[ Riley (2001) Buhler (2009) 


Vladimirov (2010)| 


Vladimirov (2011) There are four major incentives for the research 



described in this paper: 

(1) To derive the averaged equations of non-degenerated oscillatory flows. Such flows 
appear in various applications and their study represents a necessary step towards the 
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systematic studies of various classes of degenerated flows. 

(2) To make all the calculations and derivations elementary and free of any physical and 
mathematical assumptions (except the most common ones, like the existence of a differ- 
entiablc solutions). To use only Eulerian description and Eulcrian averaging operation, 
since they are the most transparent and allows to avoid additional steps and difficul- 
ties. For example, Eulerian description does not experience well-known difficulties with 
chaotic trajectories. Such a deliberately simple framework would allow to build further 
physical models based on a solid ground. 

(3) To study the reduction of the Reynolds-stress-type terms to the drift-type terms, see 



Craik fc Leibovich (1976)||Riley (200l)|| Vladimirov (20f 0)[|Vladimirov (20"TT)||Ilin fc Morgulis (20f"T) 



This rather surprising reduction was discovered by Craik & Leibovich (1976) in the small- 
amplitude theory of Langmuir circulations, but its generality and importance in many 
other areas of fluid dynamics has not yet been recognized. In particular, it is interesting 
to find what else (besides the drift velocities) can be produced by the Reynolds-stresses- 
type terms. 

(4) To understand how universal is the structure of the averaged equations. In order to 
describe universal features and variations of transformations of the Reynolds-stresses- 
type terms (or cross-correlations) we study three different problems: (i) a passive scalar 
admixture, (ii) a passive vectorial admixture (magnetic field in kinematic MHD), and (iii) 
an active vectorial admixture (vortex dynamics) . This target is methodologically impor- 
tant: after establishing the universal building blocks of the theory one can use them 



in various degenerated cases; examples of such use can be found in Vladimirov (2010) 
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For the averaging and transforming of the governing equations we employ the two- 
timing method, see e.g. |Nayfeh (1973)] |Kevorkian fc Cole (1996) | Our version of the 
method represents an elementary, systematic, and justifiable procedure introduced by 
Vladimirov (2005)[ |Yudovich (2006)[ |Vladimirov (2008)[ |Vladimirov (2010)| This proce- 
dure is complemented by consideration of the distinguished limits, that develops results 
of Vladimirov (2010) Vladimirov (2011) We also actively use the properties of com- 
mutators, which allow us to simplify analytical calculations. The employment of these 
methods and tools produces recurrent equations for both the averaged and oscillating 
parts of unknown fields. In all problems (i)-(iii) we present the calculations of three 
successive approximations (zeroth, first, and second). 

In Sect.2 the notations are briefly introduced and the list the main definitions (of the 
averaging operation, etc.) is presented. 

In Sects. 3-5 we show that in all problems (i)-(iii) the averaged equations exhibit 
the universal structure: the zeroth approximations repeat the structure of the origi- 
nal equations; in the first and second approximation the Reynolds-strcss-type terms are 
transformed into drift velocities and pseudo-diffusion, which are the same for all three 
problems; two additional terms (reminiscent of the mean-field terms in turbulence, sec 
Moffatt (1978) Moffatt (1983)) also appear for magnetic fields and vortex dynamics. In 
particular, the averaged motion of a passive scalar admixture is described only by a drift 
in the first approximation, and by a combination of a drift and pseudo-diffusion in the sec- 
ond approximation. It is remarkable that all mean-field terms (including pseudo-diffusion) 
are expressed by invariant operators (Lie-derivatives) that measure the deviation of some 
tensors from their 'frozen-in' values. 

Sect. 6 is devoted to the notion of a distinguished limit. Here we explain why for the 
non-degenerated velocity fields a slow time- variable ,s coincides with physical time t. Such 
a coincidence allows us to build a theory using only two time- variables s — t and r = at. 

Sect.l expands the area of applicability of all previously exposed results. We show 
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that for all three considered problems the high-frequency asymptotic solutions and the 
small-amplitude asymptotic solutions can be mathematically identical. This property is 
valid for s-indepcndent velocity fields. 

Sect.8 contains physical analysis of the nature of pseudo-diffusion. Our consideration 
suggests that this nature is purely kinematic, it is not related to physical diffusion. 

Sect.9 (Discussion) is devoted to the restrictions and limitations of the introduced 
transformations of the governing equations, as well as to their possible generalizations. 



2. Used Functions and Notations 

The variables x = (x\, X2, £3), t, s, and r in the text below serve as dimensionless carte- 
sian coordinates, physical time, slow time, and fast time. The used definitions, notations, 
and properties are itemized with bullets (•): 

• The class H of hat-functions is defined as 

feW: f(x,s,T) = f(x,s,T + 2n) (2.1) 

where the r-dependence is always 27r-periodic; the dependencies on x and s are not 
specified. 

• The subscripts t, r, and s denote the related partial derivatives. 

• For an arbitrary / £ H the averaging operation is 

(f) = 7T f(x,8,T)*r, Vr (2.2) 

where during the r-integration s = const and (/) does not depend on tq. 

• The class T of tilde-functions is such that 

/£T: f(x,s,T)=f(x,s,T + 2ir), with (/) = 0, (2.3) 

Tilde-functions are also called purely oscillating functions; they represent a special case 
of hat-functions with the zero average. 

• The class B of bar-functions is defined as 

/£»: / T =0, f(x,s) = (f(x,s)) (2.4) 

• Any H-function can be uniquely separated into its bar- and tilde- parts with the use 
of dO): 



/ = / + / (2.5) 

• The T '-integration (or the tilde-integration): for a given / we introduce a new function 
f T called the T- integral of /: 

f = J /(», s, P)dp~^ jj (£ f(x, s, p) dp) dp (2.6) 

• The unique solution of a PDE inside the tilde-class 

f T = df/dT = => / = (2.7) 



follows from (|2.6|) . 

• The r-derivative of a tilde-function always represents a tilde-function. However the 
r-integration of a tilde-function can produce a hat-function. An example: let us take 
4> = o sinr where 4> be an arbitrary bar- function: one can see that (<f>) = 0, however 



4 V. Vladimirov 

(J„ r (j>(x, s, p)dp) = <f) ^ 0, unless <p Q = 0. Formula (|2.6[) keeps the result of integration 
inside the T-class. 

• The T-intcgration is inverse to the r-diffcrentiation (/ T ) T = (f T ) T = f', the proof is 
omitted. 

• The product of two tilde-functions / and g represents a hat-function: fg = F, say. 
Separating its tilde-part we write 

F = F-(F)=fg-(fg)={fg} (2.8) 

where the notation {■} is introduced to avoid two levels of tildes. 

• A dimensionlcss function / = f(x,s,r) belongs to the class O(l) 

/eO(i) (2.9) 

if / = 0(1) and all partial x-, s-, and r-derivatives of / (required for our consideration), 
are also 0(1). 

• All large or small parameters in this paper are represented by various degrees of a large 
parameter a (which represents dimensionless frequency) only; large or small parameters 
appear as explicit multipliers in all formulae, while all functions always belong to O(l)- 
class. 

• The commutator of two vector-fields / and g is 

[f,g] = (g-V)f-(f-V)g, (2.10) 
It is antisymmetric and satisfies Jacobi's identity for any vector-fields /, g, and h: 

[f,9} = -[g,f], [f,[9,h}] + [h,[f,g]} + [g,[h,f}} = (2.11) 

• From (|2.5I) and (|2 .4[) one can see that 

fr = fr, (fr) = (fr) = (2.12) 

• As the average operation (|2.2[) is proportional to the integration over r, the integration 
by parts yields 

{frgh) + (fdrh) + (Jghr) = (2.14) 

(fg T ) = {(f T 9 T )r) - (Fg) = -(f T g) = -(.f T g) (2.15) 

<[/,&]} = -<[/r,fl]> = -<[/r.fl]>. (\f,r\) = -HTM = -HTM (2.16) 

Many of the above definitions and terms are different from the ones used in various 
branches of physics and fluid dynamics. We have introduced our own terms in order to 
avoid ambiguities. 



3. Scalar Admixture 

3.1. Two-timing formulation 

The equation for a passive scalar admixture c transported by a given velocity field u* in 
an incompressible fluid is 

Be 

— +m*-V*c = 0, V*-m* = (3.1) 

at* 

where asterisks mark dimensional variables, i*-time, x* = (x*, Xj, Xg)-cartesian coor- 
dinates, V* = (d/dx* ,8/8x2,9/8x3). This paper id focused on the transformations of 
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the governing equations, therefore the geometry a flow domain and particular boundary 
conditions can be introduced at later stages of the research (see Discussion). We accept 
that the considered class of oscillatory flows u* = u*(x*,t*;a*) possesses independent 
characteristic scales of velocity U, length L, time T, and frequency a* 

U, L, T, a*- Ti=L/U, T = T t (3.2) 

where Tj is a dependent (intrinsic) time-scale. It will be shown by the distinguished limit 
consideration that for non-degenerated flows the scales T and are of similar order 
(one can accept that T = kTi with a constant k ~ 1; however without any restriction 
of generality this constant can be taken as k — 1). Then dimcnsionless (not asteriatcd) 
variables and frequency are 

x = x*/L, t*=t/T, u = u*/U, <t = o-*T; 1/a < 1 (3.3) 

where 1/er is our basic small parameter. Due to the linearity of (|3.1[) we consider the 
scalar field c being dimensionless. The dimensionless version of (|3.1|) is 

c t +w-Vc = (3.4) 

The given oscillating velocity u is taken as a hat-function (|2.1[) 

u = u(x,t,at) or u = u(x,s,t) with r = at, s = t (3-5) 

where s and r are two mutually dependent time-variables (s is slow time and r is /asi 
time). We assume that the solution of (|3.4[) . (|3.5p also represents a hat- function 

c=c(a;, s,t) (3.6) 

In order to justify this assumption one needs to build such an unique solution for arbitrary 
initial data. The chain rule transforms (|3.4[) to 

d 1 d \ 1 

^- + ~T- )c+-u-Vc = (3.7) 
or a as J a 



Eqn. (|3.7[) contains the only small parameter 

1 1 



Tct* a 



(3.8) 



Here one must make an auxiliary (but technically essential) step: after the use of the chain 
rule p.7p the variables s and r arc considered to be (temporarily) mutually independent: 

t and s are independent variables (temporarily) (3-9) 

From the mathematical viewpoint the rise of the number of independent variables in a 
PDE represents a very radical step, which leads to an entirely new PDE. This step has 
to be justified a posteriori by showing that the error of the obtained solution (rewritten 
back to the original variable t and substituted into the original equation (|3.4[l ). is small. 
Hence there is the equation which is going to be considered in detail 

c r + e(u ■ V)c + ec s = 0, divw = e = l/a^0 (3.10) 

3.2. Results for the successive approximations of scalar field 
We look for the solutions of (|3.10|) in the form of regular series 

oo 

(c,«) = ^e fe (^,u fc ); c kl u k eHnO(l), k = 0,1,2,... (3.11) 

k=Q 
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where all Uk are given. The main approximation of the general velocity field u is not 
degenerated: 

Uq = uq + uq, with Ho ^ and uq ^ (3.12) 

The substitution of (|3.11[) into (|3. 10[) produces the equations for the first three successive 
approximations 

?0r = (3.13) 
ci T +u -Vc + c 0s = (3.14) 
c 2r + ui ■ Vcb + So • Vci + cis = (3.15) 
C3r + "2 ■ Vcb + So • Vc 2 + Mi ■ Vci + c 2s = (3.16) 

c fc = Ck(x,t) + c k (x,t,T), c k eBnO(i), c k e TnO(i), fc = 0,1,2,3 

The detailed solution/transformation of (|3.13p - (|3.16j) is given in Appendix A. Here we 
formulate the result. The truncated general solution is 

= c + eci + £ 2 ?2 + e 3 c 3 (3.17) 

Its bar-parts Cfc satisfy the equations 

c Os +M O -Vc o = (3.18) 
c ls + uq ■ Vci + (mi + V ) • Vc = (3.19) 
c 2s + M • Vc 2 + (mi + V ) ■ Vci + (3.20) 

(^ + Fi 0+ Fi + F x ).v,i = A(^g) 

F = i<[a ,€]>, = i<[[So,€U]), (3.2i) 

Vio = ([«!,€]), V x = i<K,£C]), S, fc ^£<6&>, (3-22) 
where the operator £ (Lie-derivative) is such that its action on £ is defined as 

Z£=L + [l,uo] (3.23) 

and its action on any (constructed from £) tensorial field fik... is such that ££ = implies 
£/ifc... = 0. In particular, 

= (ds + mo • V) - ^ fe£ m ) - (£ fc £ m ) (3.24) 

OXfn OX m 

Formulae (|3.23[) and (|3.24[) are also known as 'frozen-in' operators for a vectorial field 
and tensorial field. After (|3.18[) - (|3.20[) are solved, the tilde-parts Cfc of (|4. 1 1 1) are given by 
the recurrent expressions 

c a = 0, (3.25) 

ci = -i-Vc , (3.26) 

c 2 = -u\ • Vc - mo • V2I - \ ■ Vci - {m • VHi} T - Z[ s , (3.27) 

C3 T = — M 2 • Vco — Uq ■ Vc 2 — Mo • Vc^ — u[ ■ Vci — Ml • Vc[ — (3.28) 
{M -Vc 2 } T -{Si-VHi} r -c^ s 

Three averaged equations (|3.18j) - (|3.20j) can be written as a single advection-pseudo- 
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diffusion equation (with the error 0(e 3 )) 

d ( dc\ 

(d s +v-V)c= — U lk — j (3.29) 

c = = c Q + £ci + e 2 c 2 , v = + V l2] (3.30) 

u [2] = u a + eux + e 2 u 2 , V [2] = eV + e 2 {V 10 + Vi + V x ) (3.31) 

Kik = xfk = ? 2 Xik (3-32) 



where all the coefficients are given in p.21j) . (|3.22p . Eqn. (|3.29p shows that (with a given 
precision 0(e 2 )) the motion of c represents an advection with velocity 17 and pseudo- 
diffusion with matrix coefficient Kik = £ 2 Xik- 
One can write the truncated solution 

=c Q + c 1 +c 1 +c 2 +c 2 + c 3 + c 3 (3.33) 

where all the terms except C3 arc determined by the above equations. After rewriting 
(|3.33|) back to the original variable t (s = t,T = at) and its substitution into (|3.4j) one 
can show that the error term in the RHS of the equation is 0(e 3 ). This estimation might 
be considered as mathematical justification of the used procedure, including its most 
sensitive part (|3.9p . 

4. Kinematic MHD- equations 

4.1. Two-timing formulation of the kinematic MHD 

The equation for a magnetic field h to be 'frozen' into a given oscillating incompressible 
velocity field u*(x*,t*;a*) is 

— + [h,u*}*=0 1 V*-m*=0, V*-h = (4.1) 

where all the notations are the same as in Sect.3.1; [-,-]* stands for dimensional com- 
mutator (|2.10[) , and h is taken dimensionless due to linearity. The analysis of dimensions 
is the same as in Sect.3.1; the dimensionless version of (|4.ip is 

h t + [h, u] = 0, V • u = 0, V • h = (4.2) 

The given oscillating velocity is taken as a hat-function (|2.ip 

u = u(x, s, t); with r = at, s = t (4-3) 

and the solution of (|4.2p . (|4.3p is also taken as a hat-function 

h = h(x 1 s 1 r) (4.4) 

Then the following equation has to be studied 

h T + e[h, u] + eh s = 0, e^O (4.5) 

where the variables r and s are again temporarily independent (|3.9I) . 

4.2. Results for the successive approximations of a magnetic field 
We look for the solution of (|4.5p in the form of regular series 

00 

(h,u) = Y, £k (hk,u k ); ^,2 fe 6lnO(l), k = 0,1,2,... (4.6) 

k=0 
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where all u k are given, and (as in Sect.3.2) we consider a prescribed non-degenerated 
velocity field: Uq ^ and Uq ^ 0. The substitution of (|4.6[) into (|4.5p produces the 
equations for the first three successive approximations 

h 0T = (4.7) 
h lT + [h ,u a }+h s=0 (4.8) 

ft2 T + hl s + [ft ,«l] + [^l,«o]=0 (4.9) 

h 3r + h 2s + [h , u 2 ] + [h 2 ,u ] + [hi,ui] = (4.10) 
h k = h k (x,t) + h k (x,t,r), h k eBnO(l), h k eTnO(l), k = 0,1,2,3 
The detailed solution/transformation of (|4.7p - (|4.10p is given in Appendix B. Here we 

-[3] 

formulate the result. The truncated general solution h is 

h = h + ehi + e 2 h 2 + e 3 h 3 (4.11) 
The bar-parts h k satisfy the equations 

h 0s + [h ,u } = (4.12) 

h ls + \h o ,u 1 + V } + \h 1 ,u o ) = (4.13) 

h 2s + [h 2 ,u ] + [Tii, (ui + Vo)] + [ho, (U2 + V w + Vi+V x )] = (4.14) 
d (_ dh 0i \ -r dh M _ T 

were the notations are the same as in (|3.21|) . (|3.22|) and 

b ikl a ik =£U l g^-——) (4.15) 

The operator £ is Lie-derivative, see ([3~23]) , ([334]) . After (|432> (|4T4|) are solved, the 
tilde-parts h k of f|4. 1 1[) are given by the recurrent expressions 

h Q = 0, (4.16) 

fci = [55,ho] (4.17) 

h 2 = -h ls - [h ,ul] - [h 1 ,u ] - [hi,ul] - {[hi,u ]} T (4-18) 

h 3 = -h 2a - [h ,u 2 ] - [hl,u ] ~ [h 2 ,ul] - (4.19) 

[hi, «i] - [hi, Si] - {{h 2 ,Uo}} T - {[h 1 ,tti]} r 

Three averaged equations (|4.12p - (|4.14p can be written as a single equation that combines 
advection, pseudo-diffusion, and other mean-field- type terms (with the error 0(e 3 )) 

h = h m = ho + ehi + e 2 h 2 (4.21) 

The notation for v, vP\ and Ki k arc the same as in Q3.30p . while B ik \ = e 2 bi k i 

and Ai k = e 2 ai k . Eqn. (|4.20p shows that the evolution of h represents (with the error 
0(e 3 )) an advection and stretching by velocity v, pseudo-diffusion with matrix coefficient 
Ki k = s 2 Xi k , and other mean- field deformations. The terms with B^i and Ai k are essen- 
tially new in comparison with (|3.29j) . The construction and justification of the truncated 
solution can be performed in a similar to Sect.?) way. 
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5. Euler's equations 

5.1. Two-timing problem of vorticity dynamics 

The governing equation for dynamics of an inviscid incompressible fluid with velocity 
field u* and vorticity w* is taken in the vorticity form 

^- + [u>*,u*]* =0, V*-u*=Q (5.1) 
u*=V*xu* (5.2) 



where all the notations are the same as in Sects. 3.1 and 4.1. The equations (|5.ip arc 
mathematically the same as (|4.1j) in MHD-kincmatics. However, they are complemented 
by an additional constraint (|5.2[) . As a result the velocity field represents an unknown 
function and can not be considered as the prescribed one. It is remarkable, that this very 
sound difference does not have any impact on all the results of Sect A and the derivations 
of Appendix B. The sufficient (for the derivation of the averaged equations for vorticity) 
operations are: we have to replace ft, by a; and write the constraint (|5.2[) additionally to 
all equations. Of course, the same averaged equations can be obtained independently. 
The system of dimcnsionlcss equations under consideration can be written as: 

T + e[w, u] + eui s = 0, w = V x u, div u = 0; e —> (5-3) 

We are looking for the solution of these equations in the form of regular series 

oo 

(0,2) = ^e fe (w fe ,fi fc ), Q k = Vxu k - w fe , u k e H n O(l), k = 0, 1, 2, . . .(5.4) 

fc=0 

where all u k represent unknown functions. Here we formulate the result. The truncated 
general solution O' 3 ' is 

= O + eui + e 2 Q> 2 + e 3 ll> 3 (5.5) 

The bar-parts ZJ k satisfy the equations 

u>0s + [uj ,u ] = (5.6) 

wi s + [5Jo,«i + V ] + [wi,u ] = (5.7) 

uJ 2s + p2,u ] + [ZS t , (ui +V )} + po, (u 2 + V w + Vi +V X )} = (5.8) 
d (_ duJ 0l \ - duJoi _ _ 

All these equations are complemented by the constraints uJ k = Vxu k for k = 0, 1, 2. After 
(|5.6j) - (|5.8j) are solved, the tilde-parts u> k of (|4.11j) are given by the recurrent expressions 

Oo = 0, (5.9) 
wi = [ul,U Q ] (5.10) 
2 = -0[ s - [wo, i*i ] ~ [w[,m ] - [wi,Mq] - {[wi,m ]} t (5-11) 

0.3 = -W2 S - [u>0,ul] - [W2,M0] - [w"2,«5] - (5.12) 

[01, Iti ] - [wi,u[] - {[O 2 ,w ]} T - {[Oi,mi]} t 



All these equations are complemented by the constraints u) k = V x u k for k = 0, 1,2. 
Three averaged equations (|5.6j) - (|5.8|) can be written as a single equation that combines 
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advection, pseudo-diffusion, and other mean- field- type terms (with the error 0(e 3 )) 

— . r 1 d (- ^oA t dujm - - „- 1Q , 

w s + [w,«J = — I Xm-q^t J ~ biM ~Q^ _ aifeWfe (5.13) 

v = u Q +e{u! +V )+s 2 (u 2 + V w + V 1 + V X ) (5.14) 

^fe=x[fe=e 2 Xifc (5-15) 
uJ = u' 2 ' = uJ + ewx + s 2 uj 2 (5.16) 

Eqn. (|5 . 1 3[) shows that the evolution of tJ represents (with the error 0(e 3 )) an advection 
and stretching by velocity v, pseudo-diffusion with matrix coefficient ~Kik = £ 2 Xiki an d 
other mean-field deformations. It is very important that uJ ^ V x U. 
One can now write the truncated solution 

lD' 3 ' = ZU + u>i + u>\ + u>2 + U2 + ^3 + w 3 (5-17) 

where all the terms except SJ3 are determined by the above equations. Its mathematical 
justification is similar to the one presented in Sect.3. 

The difference between the results for vortex dynamics and the results for a magnetic 
field is essential: in the former case velocity is prescribed, while (|5.6|) - (|5.8j) represent 
the system of equations for unknown velocity and vorticity. The only exception is the 
main term uo of oscillatory velocity, which is potential (due to (|5.9[) ). and remains to 
be prescribed. There are several ways to prescribe it in different problems: (i) it can be 
forced by oscillating boundary conditions; (ii) it can appear as self-oscillations; (iii) it 
can be maintained by an external oscillating force; or (iv) it can appear in full viscous 
theory from the procedure of matching an outer flow with an oscillatory boundary layer 
solution. Hence two terms of drift velocity Vo and Vi (|3.2ip represent the functions 
that are 'external' for the averaged equations (|5.6[) - (|5.8|) . while all other terms of drift 
velocity and all mean-fields (including pseudo-diffusion) are explicitly defined by both 
Uq and unknown solutions of (15.61) . 



6. Distinguished limit 

We formulate the idea of a distinguish limit using the problem of Sect.3 (scalar trans- 
port) for a special (s-independent) non-degenerated velocity field 

u = u(x,t), 2eO(l)nH, T = at, s = l/a (6.1) 

Let us select the class of trial slow time- variables as s — a a t, where a = const (it should 
be a < 1 since s-variable must be 'slow' in comparison with r = at). Then the two-timing 
form of eqn. (|3.4j) is 

( ^ + e 1 - a 4-)c + eu-Vd=0 1 c=c(x,s,r) (6.2) 

\OT OS J 

where r and s are independent variables. Surprisingly, the slow time-scale (and the value 
of a) is uniquely determined by the structure of the equation. It is determined as the 
distinguished limit, which includes the following properties: (i) the solution for a = ad 
is given by a valid asymptotic procedure; (ii) all solutions for ad < a < 1 contain terms 
secular in s, and (iii) for any a < ad the system of equations for successive approximations 
contains internal contradictions and it is unsolvable (unless it degenerates). One can show 
that in all three problems, considered in Sects. 3-5, ad — 0. We have fully studied the first 
three approximations for all three problems, therefore it can be accepted that the validity 
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of the procedures with ad = (item (i) above) has been established. The properties (ii) 
and (iii) for ay = we accept without proof. The value ad = justifies our choice T = T t 
in (jOl . 

We should emphasise that the nature of time-variables t and s is different: the 'mod- 
ulation' variable t in u(x,t,cjt) p.5[) is given by the original physical formulation of 
the problem, while slow time-variable s in the equation p.lOp uniquely appears from 
the distinguished limit consideration. All considerations of Sects. 3-5 are based on the 
coincidence (which is the same as ad = 0) of slow time-variable s with physical time t. 
This lucky coincidence allows us to return from a special formula for velocity (|6.1|) to the 
original general functional form (|3.5[) . In more general situations (when the scales s and 
t have different orders) one should consider the problems with three time- variables: s, t, 
and r. In our case s = t, therefore we do not need to use the three-timing method. 

One more example of a distinguished limit can be seen in the averaged equations 
(|3.18p - (|3.20p . Let us take a degenerated (predominantly oscillating) flow Ho = 0, uq ^ 
and look for a solution with a = 0. In this case eqns. (|3.18j) . (|3.19p give 

cos = (6.3) 
c ls = -(M 1 + y )-Vc (6.4) 

The integration of the first equation gives Co as an arbitrary function of x only; its 
substitution into the second equation shows that (in the generic case) c\ grows linearly 
with s. The roots of this trouble lie in the fact that for predominantly oscillating flows 
ad = — 1 and s = t/a, see Vladimirov (20lT)| a wrongly chosen variable s inevitably 



leads to the appearance of secular terms. 



7. Conditional isomorphism between the high-frequency asymptotic 
solutions and small-amplitude asymptotic solutions 

The area of applicability of the derived above results is essentially expanded by the 
isomorphism between the high-frequency asymptotic solutions and the small-amplitude 
asymptotic solutions. In order to clarify the existence of such isomorphism, let us rewrite 
eqn. (|3.4p for a newly defined unmodulated velocity u 

U = £ 1 Ui(x,T), Mi G 0(1) l~l H, T = cr£, £7 = 0(1) (7.1) 

c T + ei-ui ■ Vc = 0, divui^O, ei -> (7.2) 
a 

where s\ is a small amplitude of velocity u. The most important difference with the case 
of Sects. 3-5 is a = 0(1). The two-timing form of solution is 

c = c(x, si,n) where n = r; si = e^r; = + e{ (7.3) 

OT OTl OSl 

where Si is a trial slow time- variable with an indefinite constant /3 > 0. The substitution 
of dLSJ) into ([721) gi ves 

^r- + £ i7r-W £ i- s i- v?=0 ( 7 - 4 ) 

OTi osi J a 

which is our stage to declare that the variable T\ and si are temporarily independent 
(c/. with (13.910 . One can see that the mathematical form of cqn. (|7.4j) coincides with 
(|6.2j) . Therefore, according to the result of the previous section, the distinguished limit 
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corresponds to ft = 1 and the final equation to be solved is 

Ctj + £ic Sl + ei— tii • Vc = (7.5) 
cr 

The comparison between (|7.5[) and p.lOp shows the mathematical isomorphism between 
the high-frequency asymptotic problem and the small-amplitude asymptotic problem, 
where the replacements 

r -H- Ti, s -H- si, m o ui/cr, e <h> ex (7-6) 

provide the transformation of the equations (|7.5|) and (|3. 10|1 into each other. 
We need to add some further comments: 

1. The slow variable Si = Si<rt is different from t, therefore, in contrast with (|3.5j) . t 
can not included to the list of original variables of velocity (|7.1|) . 

2. The area of applicability of the introduced isomorphism is restricted by velocity 
fields in the form (|7.1|) . Physically, the expressions in the form (|7.1|) mean that we are 
allowed to consider only steady velocity fields u(x) with superimposed unmodulated 
oscillations u(x,t). Such functional form can be deliberately chosen for the first two 
problems (scalar transport and kinematic MHD). For vortex dynamics one should impose 
the same functional restrictions on Uq. 

3. The existence of similar to (|T.6[) isomorphisms should be checked separately for each 
new equation. In particular, the adding of diffusivity or viscosity can alter the appearance 
of isomorphism. 

4. One can find that the discussed isomorphism can be also established by the simul- 
taneous rescaling of time- variable r and velocity u in (|7.1[) . We think that our way is 
preferable since: (i) it preserves the original physical variables; (ii) it can be used in more 
general equations. 

5. The high-frequency problems and small-amplitude problems are physically different. 
The introduced isomorphism shows only their mathematical equivalence, which is valid 
under the formulated restrictions. 

6. The existence of isomorphism (|7.6p explains why some high-frequency solutions and 
small-amplitude solutions can coincide (or can be close to each other). For example, 



it explains why the average equation for vorticity by Craik & Leibovich (1976) is very 
similar to eqn. (|5.8p and represents a steady version of the MHD-Drift equation from 
Vladimirov (2011)] taken without magnetic field. 



8. On the nature of pseudo-diffusion (PD) 

The pseudo-diffusion (PD) term that appears in the averaged equations (|3.20p . (|3.29l) 
might be seen as a physically new result. The prefix pseudo- appears in our study due 
to three reasons: (i) in the successive approximations PD appears for the first time in 
the equation of the second approximation (not in the zero-order equation); hence, math- 
ematically PD always plays a part of the known RHS in the equation of the second 
approximation; (ii) we use only regular asymptotic expansions (|3.11l) : and (iii) PD coef- 
ficient (matrix) has a specific form of Lie-derivative (|3.20p . (|3.24p . The first two reasons 
are very common (see e.g. Magar & Pedley (2005) I, they comply with the role of PD as 
an effective physical diffusion. However, the third reason might lead us to the suggestion 
that this term gives a s-dependent kinematic deformation of the averaged field, which is 
not related to physical diffusion. We illustrate the role of this term in a simple example 
of three-dimensional translational rigid-body oscillations. 

Let an infinite fluid oscillate translationally as a rigid-body with a small displacement 
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x(s, t) = O(e), where e = I /a. Then the Eulerian coordinate of a particle is x = x + x. 
The drift velocity identically vanishes in all orders V = 0, hence for any particle x = 
const. Therefore, we accept that 5(0,0) = 0, which means that the Lagrangian (initial) 
coordinate of each particle coincides with its unperturbed position X = x. Then the 
solution is 

c(x,s,t) = c L (x -x) (8.1) 

where c L (X) is a given time-independent distribution in Lagrangian coordinates. One 
can expand both sides of ()8.1[) as 



^ t z \ ~dc L (x) e 2 _ d 2 c L (x) , . 

c Q + EC X + e 2 c 2 + ... = c L (x)- Xi ^ + —XiX k + . . . (8.2) 

where the LHS corresponds to the decomposition (|3.1ip with c„ = c„(x, s, t) = c n (x, s) + 
Cn(x, s, t), while the RHS represents the Taylor's series. The averaging (|2.2I) of (|8.2|) yields 

co + &h + e 2 c 2 + ... = c L (x) + ^(Mk) d g^ +■■■ (8.3) 

where we have changed x(s, r) to e£(s, r) p.21[) . which is valid with the given precision 
(by an elementary consideration one can derive that x = e$, — e 2 £ s + 0(e 3 )). The s- 
derivative of (|8.3p is 

-co, + s-c ls + e 2 -c 2s + ... = ^(6&>.^ + (8-4) 

where we have used c L = cq (which is the leading term of (|8.3[) 1 . One can see that 
cqn. (|8.4j) coincides with (|3.18|) - (|3.20|) taken for Ho = and for the zero drift. 

This example indicates that PD represents a purely kinematic correction which appears 
due to the averaging operation. In addition, one might put forward a conjecture that all 
mean-field terms in (|3.29[) . (|4.20j) . (|5.13|) are also purely kinematic and do not represent 
any physical processes. However, the same 'philosophical' question arises in the theory 
of turbulence: let us take an ensemble of flows, which differ from each other only by 
translations as a whole; then the averaging operation with a time-dependent distribution 
function leads to 'turbulent diffusion'. Can it be considered physically meaningful? The 
answer to this question requires further analysis. Anyway, the significance of PD (and 
other mean-field terms) for the results of this paper lies in the discovery of the universal 
mathematical structure of the terms, that appear from the straightforward calculations 
in addition to the drift-type corrections of the averaged equations. 



9. Discussion 

1. The main novel element of our consideration is a universal structure of averaged 
equations describing different processes in non-degenerated oscillatory flows. The most 
useful problem (for the analysis of this structure) is the one for a scalar admixture 
(Sect. 3), where the calculations are the simplest and at the same time sufficient enough 
to produce the analytical expressions and constructions, that appear in all other cases, in- 



cluding predominantly oscillating (degenerated) flows, see Vladimirov (2010) Vladimirov (2011) 

2. The averaged equations for a 'frozen-in' magnetic field and vorticity dynamics for- 
mally coincide; the difference between two cases is the constraint linking velocity with 
vorticity in the latter one. However, the presence of such a constraint has entirely changed 
the nature of the equations: the velocity field is prescribed in the MHD-kincmatics but 
it represents an unknown function in vortex dynamics. 
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3. For all considered non-degenerated flows a non-linear averaged equation appears 
only in the mean (zero) approximation of vortex dynamics; its form coincides with that 
of the original (exact) equations. All other equations arc linear. The most impressive 
in the analytical beauty and demanding physical explanations is an appearance of Lie- 
derivatives in the pseudo-diffusion matrix and in other mean-fields. 

4. In this paper we intensely use the notion of a drift. It is known that a drift ve- 
locity can appear from Lagrangian, Eulerian, or hybrid (Euler-Lagrange) considerations. 
In our study we use the Eulerian drift, which appears as a result of the Eulerian av- 
eraging of the governing PDEs without direct addressing the motion of particles, see 



Craik fc Leibovich (1976)l|Craik (1985j||Riley (200 1)|| Vladimirov (2010)[|Vladimirov (20lT)| 



Ilin & Morgulis (2011) More classical Lagrangian drift appears as the average motion of 



Lagrangian particles and its theory is based on the averaging of ODEs, see Stokes (1847) 
Lamb (1932) Longuet-Higgins (1953) Batchelor (1967)| the hybr id drift coincides with 



the Lagrangian one, see |Andrews fc Mclntyre (1978)[|Craik (1985)||Soward fc Roberts (20"l~0) 
Vladimirov (2010) It is known that Lagrangian and Eulerian drifts are similar, but not 



identical to each other, see Vladimirov (2010) 

5. The first two examples (the transport of a scalar and the kinematic MHD) are not 
always mutually independent: for the translationally- invariant motions eqn. (|4.1[ ) can be 



reduced to (|3.ip . see Vladimirov, Moffatt and Ilin (1996) 



6. The small parameter of our asymptotic theory is inverse dimensionless frequency 
e = 1/oj. One can see that (since the oscillatory velocity is 0(1)) the dimensionless 
characteristic oscillatory displacement of a material particle is 0(e). It is interesting, 
that for the isomorphic small-amplitude flow of Sect.7 this displacement is O(ei). 

7. For a finite and time-dependent flow domain T>(t) the definition of average (|2.2|) 
works only if x G T> at any instant. If it is not true, then the theory should include the 
'projection' of a boundary condition on an 'undisturbed' boundary. Such a consideration 
requires the smallness of the oscillatory displacements of fluid particles, which are always 
proportional to e (see item 6 above). Hence, the generalizations of all our results to finite 
domains can be undertaken without adding any new small parameters to the list (|3.8j) . 

8. Viscosity or diffusivity can be routinely added to all the above considerations as the 
RHS-tcrm «V 2 a, v\/ 2 u, etc. However, the final appearance of such terms depends on 
the order of magnitude of the new small parameter (like 1/Pe = n/aL 2 , where Pe is the 
Peclet number) in terms of s. Different relations between small parameters will produce 
different averaged equations. The most 'harmless' case of adding viscosity or diffusivity 
is to take k/<tL 2 = 0(e 2 ). In this case viscous or diffusion terms will be just added to 
all the equations of the second approximation. This assumption has been intensely used 



in biological applications, see Magar & Pedley (2005) It is interesting to note that in 



such cases dimensionless physical diffusion and pseudo-diffusion are of the same order of 
magnitude. We avoid viscosity and diffusivity in this paper in order to focus our attention 
on 'pure' kinematics and dynamics. 

9. The incorporation of density stratification and/or compressibility into all considered 
problems is straightforward. In fact, a similar theory of the predominantly oscillating 



flows was developed for a compressible fluid, see Vladimirov (2010) 



10. In this paper we consider periodic (in fast time-variable r) functions. The stud- 
ies of non-periodic in t oscillations can represent the next stages of research. For ex- 
ample, such a development has been already available for Langmuir circulations, see 
Craik fc Leibovich (1976)] |Leibovich (1983)] |Craik (1985) 



1 1 . The described transformations make partially visible the same properties of the gov 
crning equations of fluid dynamics as were studied by the hybrid Euler-Lagrange GLM 
procedure, see |Andrews fc Mclntyre (1978)[ |Buhler (200~9)j |Soward fc Roberts (2010) 
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However, the GLM-method in its original form is not adapted to the two-timing frame- 
work: it is sufficient to note that a fast-time variable and a related small parameter (the 
ratio of two time-scales) do not appear in Andrews & Mclntyre (1978) The two-timing 



adaptation of the GLM-kincmatics has been developed by Vladimirov (2010) Further 



steps in this direction arc required to introduce the two-timing dynamics. 

12. Degenerated flows with u = (predominantly oscillating flows) were consid- 
ered by [Stokes (1847)[|Craik fc Leibovich (1976)[|Riley (2001)| jMagar fc Pedley (2005) 



Buhler (2009)| |Vladimirov (2010)[ |Herreman fc Lesaffre (2011)[ |Vladimirov (2011)[ and 
many others. The results of Vladimirov (2010)} " Vladimirov (2011) show that the aver- 
aged governing equations of degenerated flows contain the same (or similar) building 
blocks/elements, as the ones presented above. The vortex dynamics of such flows is ex- 
tremely interesting theoretically, since it produces non-linear equations for averaged vor- 
ticity, that include a drift in their coefficients. A classical application of this theory is the 



describing of the Langmuir circulations, see Craik & Leibovich (1976) Leibovich (1983) 



Thorpe (2004) 



13. There are many areas where the results of Sects.3-7 could be applied, such as 
waves in a flow with strong shear, flows past translationally moving and oscillating solids, 
oscillatory flows in channels and pipes, various acoustical problems, etc. The number of 
applications can be very large. However we do not consider any applications, since: (i) 
the focus of this paper is the universal structure of the averaged governing equations; (ii) 
every application of a general theory requires its own paper. 

14. This paper is written as Part 1 of a series, where Part 2 (a similar theory for 
degenerated flows) is based on Vladimirov (2010) A number of examples illustrating the 
values of drift and pseudo-diffusion is given in Part 2, since predominantly oscillating 
flows have many well-developed applications. 

The author is grateful to Profs. A.D.D. Craik, H.K.Moffatt, T.J. Pedley, and A.M.Soward 
for helpful discussions. The author thanks the Department of Mathematics of the Uni- 
versity of York for the research-stimulating environment. 



Appendix A. Deriving of the equations for a passive admixture 

The equation of zero approximation (|3.13[) is 

c Ot =c Ot = (Al) 

Its unique solution is cq = (|2.7p . while cq remains undetermined; let us write it as 

c = 0, (A2) 

where the first equality gives (|3.25j) . The equation of the first approximation (|3.14j) (with 
the use of (|A~2)) ) is: 

Ci T + u ■ Vc + c Qs — (A 3) 

Its bar-part (|2.4[) and tilde-part (|2.3I) arc 

co s + u • Vc = (A 4) 

c lT + u ■ Vc = (A 5) 

where (|A 4[) represents the final equation for cq (|3.18[) . The T- integration (|2.6I) of (|A 5|) 
yields (|3.26|) and keeps C\ unknown 

ci = -Uq ■ Vc , ci =[T| (A 6) 
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The equation of second approximation (|3.15[) that takes into account (|A 2\ is 

c 2t + «i • Vc + mo • Vci + ci s = (A 7) 

Its bar-part and tilde-part are 

c u + ui ■ Vc + M • Vci + (S • Vci) = (A 8) 

c 2r + u\ ■ Vc + «o • Vci + m • Vci + {S • Vci} + cis = (A 9) 

where (•) stands for the average (|2.2p and {■} for the tilde-part (|2.8p . Eqn. (|A 81) can be 
transformed (with the use of (|A6j) and ([230]) - ([236]) ) into the form (f3~19l) 

c is + Ho • Vci + (mi + V ) • Vc = (A 10) 

V a ^ l -([u^ul]) (All) 

while the T-integration (JUJ) of (fA~9| gives (pHT)) 

? 2 = -m1-Vc -m -VS[-^-Vci-{m -VS' 1 } t -? : I s , c 2 =0 (A 12) 

The equation of third approximation (|3.16p that takes into account (|A 2[) is 

C3t + M 2 • Vc + Mo ' Vc 2 + Mi • Vci + C 2s = (A 13) 

Its bar-part and tildc-part are 

c 2s + m 2 • Vc + m • Vc 2 + ui ■ Vci + (So • Vc 2 ) + (Si • Vci.) = (A 14) 

C3t + S 2 • Vc + m • Vc 2 + u ■ Vc 2 + Si • Vci + Ml ■ Vci + (A 15) 
{S • Vc 2 } + {Si • Vci} + c 2s = 

We transform the 'Reynolds stress' terms in (|A 14|) with the use of (|A~6j) . (|A 12|) , ((2~T0| - 

R c = (Si • Vci) + (S ■ Vc 2 ) = i?cl + i? c2 + Rc3 + Rc4 (A 16) 

R cl = -((S • V)(SJ • V))ci = V • Vci (A 17) 

Rc2 = -((Si • V)(SJ • V))c - ((S • V)(S{ • V))c = Voi • Vc (A 18) 

R c3 = -((S • V){(S • V)c-i} T ) = ((S5 • V)(S • V)ci) = (A 19) 
-((uj • V)(u • V)(«J • V))c = Vi • Vc 

i?c4 = -((S • V)S[ S ) - ((S • V)(m ■ V)c[) = (A20) 

(( s S .v )M )^ x .v,„-A(x.£) 

Vo^<[«o,€]>, Vi^i([[S ,i],I]), £ = uS (A21) 

Vio s (!«!,€]>, V x = i([€,«]>, Xifcsi^fc), (A22) 
where the operator £ (Lie-derivative) is such that its action on £ and (£i£fc) is defined as 

Zl = l s + ll m ] (A 23) 

= (a, + mo • v) (Mk) - ^* (6U) - p^- (Mm) (a 24) 
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Hence eqn. (|A 14[) takes a form 

c 2s + mo • Vc 2 + («i + V ) • Vci + (A 25) 

(^+^o + F 1+ F x ).V-co = A(x, fe g) 

while (|A15|) can be T-intcgrated flU} 

?3 T = — xl 2 ■ Vco — ttp ' ^2 — wo • Vc£ — m[ ■ Vci — Mi • Vc[ — (A 26) 
{u a -Vc~ 2 Y -{ux-VZrY +^ s , c 3 =0 

Appendix B. Deriving of equations for kinematic MHD 

The equations of the zeroth approximation is 

h 0T = h 0r = (B 1) 

Its unique solution is ho = (|2.7[) , while ho remains undetermined; let us write it as 

h = 0, h =|7] (B2) 
The equation of the first approximation of (|4.5[) (with the use of (|B 2[0 is: 

h lT + [h ,u Q } +h Qs = (B3) 
Its bar-part (|2.4[) and tilde-part (|2.3[) are 

7r .s + [7r ,Mo] = o (B4) 

h lT + [^o,M ] =0 (B5) 
Eqn. (|B 4p represents a final equation for h . The T- integration (|2.6j) of ()B 5[) yields 

hx = [ulM], hi={j] (B6) 
The equation of the second approximation that takes into account (|B 2p is 

h^ T + h\ a + [h ,Mi] + [hi,Mo] =0 (B7) 
Its bar-part and tilde-part are 

h ls + [So,ui] + [Hi, M ] + ([fci,5o]) =0 (B8) 

h 2r + h is + [h ,Ml] + [hi,U ] + [ft.1, M ] + { [hi , M ] } = (B 9) 

where (■) stands for the average operation (|2.2[) and {•} for the tilde-part (|2.8[) . The 
Reynolds-stress-type term in (|B 8[) can be transformed with the use of (|B 61) and (|2.10l) - 
(j2~T6) as 

([hi, Mo]) = ([[Moh ],M ]) = [h ,V ] 

Vo^<[5o,«S]) (BIO) 

Hence (|B 8[) takes a form 

h ls + [h o ,Mi+Vo] + [hi,M O ]=0 (Bll) 
while (|B1~0|) can be T-integrated (|2~6| 

h 2 = -h ls - [h ,M[] - [hi, uo] - [hi,Mo] - {[hi,u ]} r , h 2 = \T\ (B 12) 
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The equation of the third approximation that takes into account (|B 2j) is 

h 3r + h 2s + [h ,u 2 ] + [h 2 ,S ] + [hi, Si] =0 (B 13) 

Its bar-part and tildc-part are 

h 2s + [h ,u 2 ] + [h 2 ,u ] + [hi, Mi] + ([ha, So]) + ([hi, Si]) = (B 14) 
h 3r + h 2s + [h ,S 2 ] + [h 2 ,Uo] + [h 2 ,u ] + (B 15) 

[hx.ui] + [hi, Si] + {[h 2 ,S ]} + {[hi, Si]} = 
We transform the Rcynolds-stress-type terms in (|B 14[) with the use of (|B 6I) . (|B 12|) . (|2.10[1 - 

R h = ([hi, Si]) + ([h 2 , M ]) = R h i + R h2 + R h3 + R M (B 16) 

R hl = -([[hi,S5],S ]) - [hi, V ] (B17) 

R h2 = -([[h ,S5],Si]) - ([[h , SI], S ]) = [h , Voi] (B18) 

R h3 = -([{[h!,Sar,So]) - [h ,Fi] (B19) 

-<(S5 • V)(S • V)(SS • V))c = Vi • Vc (B 20) 

i? M = -([hi, So]) - ([[hi, m ], So]) = ([£hi,S5]) = (B21) 

where the additional to (|A 2ip . (|A 22[) notations are: 

Ojfez=L(&— ), <z ifc =,L<& — — — — ) (B22) 

\ CWGj / \ OXidXk 0Xk oxi J 

Hence cqn. (|A 14|) takes a form (with the i-component in the LHS) 

h 2s + [h 2 , mo] + [hi, (mi + V )] + [ho, (m 2 + Vio + Fi + F x )]| . = (B 23) 

5 /_ <9ho*\ r dh i _ - 
^( Xlk ^- - bikl ^-~ aikhk 



while (|B 15j) can be T-integrated (|2~6| 

hs = -hL - [ho.Sa] - [h 2 ,Mo] - [h 2 ,«5] - (B24) 
[h^ tZi] - [hi, Si] - {[h 2 ,M ]} T - {[hi,Si]} T , h 3 =|T] 
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